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Abstract. The state of the matter that is obscuring a small circumnuclear region in active galactic nuclei can 
be probed by observations of its broad emission lines. Infrared lines are particularly useful since they penetrate 
significant columns of obscuring matter, the properties of which can be constrained by comparing infrared and 
X-ray obscuration. We report on new 4/j.m spectroscopy with ISAAC at the ESO VLT of a sample of 12 Seyfert 2 
galaxies, probing for broad components to the Brackett a 4.05/im hydrogen recombination line. Broad components 
are observed in 3 to 4 objects. All objects with a broad component exhibit relatively low X-ray obscuring columns, 
and our results are consistent with a Galactic ratio of 4/im obscuration to the BLR and X-ray column. In 
combination with observations of a non-Galactic ratio of visual obscuration of BLRs and X-ray obscuring column 
in Seyferts, and interpreted in a unified AGN scheme, this result can be reconciled with two interpretations. Either 
the properties of dust near the AGN are modified towards larger grains, for example through coagulation, in a 
way that significantly flattens the optical/IR extinction curve, or the ratio of dust obscuration to X-ray column 
varies for different viewing angles with respect to the axis of symmetry of the putative torus. Our spectra also 
provide a survey of emission in the [Si IX] 3.94pim coronal line, finding variation by an order of magnitude in its 
ratio to Bra. The first extragalactic detection of the [Ca VII] 4.09/im and [Ca V] 4.16/im coronal lines is reported 
in the spectrum of the Circinus galaxy. 
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1. Introduction 

Unified scenarios have been highly successful in explain- 
ing several aspects of the AGN phenomenon, by assuming 
that different manifestations of the AGN phenomenon cor- 
respond to similar objects viewed from different directions. 
The detection in polarized light of broad emis sion lines in 
Seyfe rt 2 galaxies (Antonucci & Miller 1985 , Antonucci 
1993 ) has been central to the development of these sce- 
narios. Subsequently, spectropolarimetry has become the 
prime tool for detecting hidden broad lines in larger sam- 
ples (e.g., Miller & Goodrich |i"990[ Tran et al. fL99§ Young 
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et al. |1996|; Heisler et al. |l997t Moran et al. |2000 
et al 



2001 



Lumsden 



Tran 2001). X-ray spectroscopy has been the 
second key observation, finding Seyfert 2s on average much 
more highly obscured than Seyfert Is and quantitatively 
establishing the absorbing column densities in neutral and 
'warm', i.e. highly ionized material (e.g., Turner et al. 
1997; Bassani et al. 1999). Still, relatively little is known 



about the actual distribution and physical state of the 
material obscuring the central engine of Seyfert 2 galaxies 
from our view. Is it in the form of a compact parsec scale 



torus (Krolik & Begelman 1986)? Or does obscuration on 
scales of tens or hundreds of parsec play a significant role 
(e.g. Maiolino & Rieke |1995|, Malkan et al. 1998)? Mass 



arguments suggest that at least the Compton-thick ab- 
sorbers are on scales of tens of parsecs or less (Risaliti et 
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al. 1999), but lower column components are more difficult 
to constrain. Do outflows contribute to forming the obscu- 
ration (e.g. Konigl & Kartje |L994t Elvis |2000| )? How does 
the state of the obscuring matter differ from the interstel- 
lar medium in a normal galaxy, given the extreme condi- 
tions close to a powerful AGN? High obscuring columns 
can be explained by different scenarios, and the warm dust 
emission seen in the mid-infrared cannot uniquely distin- 
guish between compact and more extended configurations 
either (e.g., Pier & Krolik |1992|; Efstathiou et al. |1995; 



Granato et al. 1997) 



One way to address some of these issues is to study 
the obscuration of the central engine in wavelength ranges 
that are partially transparent at the column densities of 
interest, and compare the results. Of particular value are 
X-rays and the infrared range, but care has to be taken to 
compare obscuration towards similar regions. Narrow Line 
Region (NLR) emission and mid-infrared dust emission 
probe regions that are much larger than the central X-ray 
source. Therefore, their obscuration may involve different 
foreground material, and cannot be compared directly to 
X-ray results. In contrast, reverberation mapping results 
show the Broad Line Region (BLR) to be well below a 
parsec m size (Netzer |1990[ ), allowing a meaningful com- 
parison of BLR and X-ray obscuration, by material that 
might be found in both parsec-scale and larger regions. 
Our goal is to constrain infrared obscuration, and thus 
the state of the obscuring matter by the detection or non- 
detection of infrared BLRs in Seyfert 2s of various X-ray 



obscuring columns. Recently, Maiolino et al. (2001a) have 
compared visual and X-ray obscuration in Seyferts, finding 
large deviations from standard Galactic values. In many 
objects, the ratio of reddening £b-v and X-ray column 
Nh is low by about an order of magnitude compared to 
the Galactic value. These results provide additional moti- 
vation to explore the relation between infrared and X-ray 
obscuration, and compare it with that in the visual. 

Practical considerations drive the choice of the optimal 
transition for Broad Line Region searches in the infrared. 
Standard interstellar extinction laws have a minimum in 
the 3-8pm range, then rise through the silicate features 
and drop again steeply beyond 30pm. This tends to ar- 
gue in favour of the longest wavelength infrared observa- 
tions. However, the flux of the strongest (a) recombina- 
tion lines drops approximately with the second power of 
wavelength, while dust continuum increases making their 
detection increasingly difficult. On the basis of such rea- 
soning and of ISO spectroscopy of Brac kett /3, Brackett a, 
and Pfund a in NGC 1068, Lutz et al. Q2000a| ) concluded 
that Brackett a with its fairly high line to continuum ra- 
tio and low obscuration is the most promising line for in- 
frared BLR searches. Using sensitive instruments such as 
ISAAC on 8m telescopes like the VLT, it is now possible 
to perform such observations and take a next step beyond 
earlier infrared searches for BLRs, which either focussed 
on lines that have a good line to continuum but still rela- 
tively high obscuration (Pa/3 1.28/im) or are suffering less 
extinction but are difficult to measure because of a low 



line-to-continuum ratio (Br7 2.17pm). Searches in these 
lines detected several broad components, but cannot probe 
beyond equivalent visual obscurations of 10-20 mag (e.g 
Rix et al 



199C; Blanco et al. 



1990 



Goodrich et al. 



1994 



Ruiz et al. |1994| Veilleux et al. |1997| ; Gilli et al. |2000D . 

This paper is organized as follows. Sect, ^discusses the 
observations and data analysis, Sect. ||tlie result of a coro- 
nal line survey obtained from our data, Sect. || presents 
the results of the Brackett a spectroscopy including line 
decompositions. We discuss the results and compare to 
X-ray and optical work in Sect. [5] and conclude in Sect. [| 

2. Observations and data analysis 

To identify suitable bright Seyfert 2 galaxies we have 
used the Bassani et al. (1999) sample with well-known 
X-ray obscuring columns. We have selected the objects 
that are brightest in extinction-corrected [O III] A5007 
(>10 _12 ergs _1 cm -2 ), observable from Paranal (declina- 
tion <20deg), and are at z<0.015 to keep Brackett a in 
the observable part of the atmospheric L band. Of the 
23 objects in this parent sample, 12 were observed fol- 
lowing right ascension and weather constraints during a 
run in March 9-11, 2001. Both the parent sample and the 
observed targets cover a wide range of X-ray obscuring 
columns. NGC 1068 was part of the parent sample and is 
included in our discussion, using the existing ISO data of 
Lutz et al. ( |2000a| ; |2000bD observed with a large 14" x 20" 
aperture. Due to the prominence of the central L band 



peak in NGC 1068 (Alonso-Herrero et al. |1998[ ), the con- 
tinuum of these data is still dominated by the central com- 
pact peak and can be analysed together with the smaller 
aperture ISAAC sample. For part of our sample, spec- 
tropolarimetric observations are available in the literature, 
directly confirming the intrinsic presence of a Broad Line 
Region. Table [j] lists this information, together with some 
basic source properties and total line fluxes measured from 
the spectra. 

We used the ISAAC-LW medium resolution spec- 
troscopy mode to cover a range of 3.93 to 4.17pm at spec- 
tral resolving power «2500. This range was chosen because 
the rapid drop of atmospheric transmission makes longer 
wavelengths useless, and to ensure coverage of the [Si IX] 
3.94pm coronal line in all spectra. While we present below 
results for this and other coronal lines obtained from our 
spectra, the primary motivation for including this line was 
to provide a reference, observed with the same instrument 
setting, for the line widths in the Narrow / Coronal Line 
Region. These can be used in cases where there is ambi- 
guity between a true BLR and relatively broad NLR com- 
ponents. The bright and compact nuclei were acquired in 
the K band and centered in the 1" slit which was oriented 
north-south. Observations were done in the chopping and 
nodding along the slit scheme suited for thermal infrared 
observations. Integration times (excluding overheads) var- 
ied between 37 and 62 minutes per source. Our strategy 
was to integrate to good signal-to-noise in the narrow com- 
ponent of Bra. Then, nondctection of a broad component 
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Fig. 1. ISAAC 4/Ltm spectra of 12 Seyfert 2 galaxies. Noise increases towards the long wavelength end because of 
increasing atmospheric opacity. 



is meaningful, since in typical Seyfert Is the broad line 
flux is several times brighter than the narrow line flux. In 
the sample of Stirpe (199C), for example, the median flux 
ratio of broad and narrow H/3 is 20 and the lowest ratio 
7, for those 14 objects that do not exceed a luminosity of 
My = -22 mag. 

Data reduction followed standard procedures in 
Eclipse and IRAF. The wavelength calibration is based 
on arc spectra and on a vacuum scale. For flux calibration 
and correction for atmospheric absorptions, we observed 
both early type and G type stars. For our particular case 
of spectra near Brackett a with sometimes good contin- 



uum S/N, we found early type standards little suited be- 
cause of difficulties to correct for their intrinsic hydrogen 
and helium absorption lines, and in some cases emission 
lines. More satisfactory correction of atmospheric absorp- 
tion was achieved using early G star spectra which had 
first been corrected for their significant intrinsic spectral 
structure using the high resolution solar spectrum pro- 
vided at the ESO ISAAC web pages. This solar spectrum 
was convolved to the appropriate resolution and slightly 
shifted and scaled, to give optimum cancellation of G star 
features when dividing a G star and solar spectrum. Since 
we are interested mostly in the inner region showing dust 
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Fig. 2. Spectra sorted in order of increasing X-ray obscuring column. To facilitate comparison, the spectra have been 
shifted to a rest wavelength scale, continuum subtracted and scaled to the same Bra peak flux. 



and BLR rather than NLR emission we used optimum spectra. Only Circinus showed evidence for faint extended 

extraction of the spectra to get the best signal-to-noise. Bra emission peaking wlO-15" from the nucleus, which we 

Effects on the NLR fluxes which in principle occur were did not include in the extracted spectrum. Figures |l]and | 

verified to be small by inspection of the 2-D spectra and show the spectra, first as observed and then sorted by the 

comparison of optimum extracted with directly extracted X-ray obscuring column and scaled in a way facilitating 
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Fig. 3. Histogram for the ratio of [Si IX]3.94/xm and the 
narrow component of Bra in 13 Seyfert galaxies 

comparison. Total fluxes of the coronal lines and of Bra 
(from direct integration of the line profiles) are listed in 
Table |l|. Due to the narrow slit and optimum extraction 
scheme, they correspond to a small ssl" aperture. 

3. A survey of coronal lines 

The [Si IX] 3.94/im line was first observed in an extra- 
galactic object by Oliva et al. ( |1994| ) in the spectrum of 
the Circinus galaxy. It is one of the highest excitation 
coronal lines observed in the spectra of Seyfert galaxies, 
with a lower ionisation potential of 303 eV. It is thus an 
interesting diagnostic of the coronal line region of AGN, 
the density of which can be constrained to be less than 
about 10 6 cm~ 3 using the d ensity- sensitive ratio [Si IX] 
2.58/3.94/im (e.g. Lutz et al. |2000b| ). Our good S/N 4/zm 
spectra of 12 Seyferts significantly widen the database of 
observations of this line. Figure || shows the distribution of 
the ratio of [Si IX] 3.94/im and the narrow (NLR) compo- 
nent of Bra, where we have again added NGC 1068 (Lutz 
et al. [2000b| ) to our sources. While the [Si IX] line is de- 
tected in 12 of 13 objects, the observed excitations dif- 
fer strongly, with the ratio of [SilX] and Bra varying by 
about an order of magnitude, similar to variations in the 
strength of optical coronal lines in Seyferts (e.g., Penston 
et al. |1984| , Erkens et al. |1997| ). 

Factors that might bias the observed ratio with respect 
to the intrinsic NLR one need consideration. Firstly, Bra 
may be contaminated by starburst emission. This effect 
is likely to be small in our small aperture. This can be 
tested on the basis of the line profiles. On the high excita- 
tion end, for Circinus both Bra and the coronal lines have 
similar profiles suggesting a common origin. Any starburst 
contribution to the observed Bra would make the intrin- 
sic NLR excitation even higher. On the low excitation end, 
both [Si IX] and Bra in Mkn 1210 have similar profiles and 
a high FWHM typical for a NLR, suggesting that both 
lines originate in the same region and that the low excita- 
tion is real. Another low excitation object, NGC 5643, has 



similar profiles for both lines but a small linewidth where 
a starburst contribution is difficult to discriminate kine- 
matically. Bra cores that are narrower than for the cor- 
responding [SilX] line appear present in NGC 3081 and 
NGC 4507. These may indicate starburst contamination 
but could also partly reflect a NLR where lower and higher 
excitation species show different dynamics, as observed 
in some Seyferts (e.g., Appenzeller & Ostreicher 1988). 
Secondly, underestimating the BLR contribution in a com- 
posite broad/narrow Bra profile leads to an underestimate 
of the [SilX] to narrow Bra ratio. In light of the Bra pro- 
files and their decomposition discussed in Sect. [|, such 
decomposition uncertainties might affect NGC 2110 and 
NGC 5506, which indeed are near the low end of the ob- 
served [SilX] to narrow Bra range. Finally, excitation gra- 
dients in the NLR may bias small aperture measurements 
like ours towards higher excitation. For the extensively ob- 
served Circinus galaxy, comparison to larger aperture and 
integral field measurements of coronal lines (Oliva et al. 



1994, Maiolino et al. 1998) suggests the possible effect to 
be less than a factor 2. 

We conclude that these effects do not dominate the 
[Si IX] /Bra spread of about an order of magnitude among 
the objects of our sample. We also computed simple 
CLOUDY photoionization models for the 'table power 
law' continuum implemented in CLOUDY and varying 
ionization parameter, showing that already modest vari- 
ations of A(log U) ~ 0.3 near log U — —2 can change 
[Si IX] /Bra by an order of magnitude. The observed [Si IX] 
flux is hence a good qualitative indicator for the presence 
of an AGN, but a relatively poor quantitative indicator of 
the AGN luminosity, due to the considerable variations of 
observed excitation. 

Two more coronal lines are detected for the first time 
to our knowledge in the spectrum of an extragalactic 
source: [CaVII] 4.09/zm and [CaV] 4.16/umin the spectrum 
of the Circinus galaxy, at wavelengths in agreement with 
the ones inferred by Feuchtgruber et al. ( 2001 ) from obser- 
vation of the planetary nebulae NGC 6302 and NGC 7027. 
Like [SilX] and Bra, they show a distinct blue asymme- 
try, the blue wing extending out to «800km/s. A hint of 
the [CaVII] line is seen in other high excitation spectra 
(NGC 4388, MCG-5-23-16) but at low significance. The 
rest wavelength of [CaV] is covered by our spectra for 
Circinus only. Detection of these lines (lower ionisation 
potentials 109 and 67eV) confirms again the unique suit- 
ability of the Circinus galaxy for coronal line studies due 
to its brightness, high excitation, and narrow line widths. 

Observations of infrared fine-structure and coronal 
lines can also help elucidating the role of dust in form- 
ing the line asymmetries and shifts that are frequently 
observed in optical coronal lines (e.g., Erkens et al. 1997). 
For the prototypical Seyferts NGC 1068 and NGC 4151, 
these issues are d iscussed on the b asis o f ISO data by 
Lutz et al. ( |2000bj ) and Sturm et al. (|l999|) . They observe 
some of the differences between optical and infrared lines 
that are expected for a moderate amount of dust in the 
NLR, but also find remaining asymmetries and shifts in 
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Source 


cz 
km/s 


[OIII] 

10 _14 Wm- 2 


N H 

10 20 cm" 2 


BLR in polarized light? 


[Si IX] 


Bra 

1Q -20 


[CaVII] 
Wm~ 2 


[CaV] 


NGC 2992 


2311 


0.680 


69 


yes (Lumsden et al. 2002) 


873 


6650 






MCG-5-23-16 


2482 


0.409 


162 


yes (Lumsden et al. 2002) 


2150 


5010 


? 




NGC 2110 


2335 


0.321 


289 




<500 


2110 






NGC 5506 


1853 


0.600 


340 


? (Tran 2001, Young et al. 1996) 


1000: 


11700 






NGC 4507 


3538 


0.158 


2920 


yes (Moran et al. 2000) 


2380 


3650 






NGC 4388 


2524 


0.374 


4200 


yes (Young et al. 1996) 


2530 


3220 


? 




NGC 4941 


1108 


0.355 


4500 


no (Moran et al. 2000) 


455 


804 






NGC 3081 


2385 


0.215 


6600 


yes (Moran et al. 2000) 


1500 


1120 






Mkn 1210 


4046 


0.482 


> 10 4 


yes (Tran et al. 1992) 


1230 


2680 






NGC 4968 


2957 


1.116 


> 10 4 




667 


1060 






Circinus 


449 


6.970 


410 4 


yes (Oliva et al. 1998) 


37700 


13900 


4040 


4610 


NGC 5643 


1199 


0.694 


> 10 5 


no (Moran et al. 2000) 


870 


2920 






NGC 1068 


1137 


15.86 


> 10 5 


yes (Antonucci & Miller 85) 


54000 


69000 







Table 1. Source propertie s and integrated line fluxes. Heliocentric redshifts are fro 
fluxes from Bassani et al. ( 1999 ). X-ray obscuring columns Nn are from Bassani et 
et al. [1999b- 



m NED, extinction corrected [O III] 
al. ( 1999 ) except for Circinus (Matt 



Source 


FWHM (b) 
km/s 


FWHM (n) 
km/s 


Av (n) 
km/s 


Flux (b) 

1Q -20 


Flux (n) 
Wm" 2 


FHWM [Si IX] 
km/s 


Av [Si IX] 

km/s 


NGC 2992 


1850 


230 


-23 


5370 


1300 


240 


-14 


MCG-5-23-16 


1450 


230 


+84 


3930 


1110 


190 


+77 


NGC 2110 




618 


-20 


<3100 








NGC 5506 


1217 


fixed: 460 


+39 


7150 


5260 


640: 


+ 150: 


NGC 4507 




540 


-13 


<4000 




520 


+74 


NGC 4388 




260 


-22 


<3100 




190 


-5 


NGC 4941 




260 


+0 


<1600 




130 


+16 


NGC 3081 




260 


-16 


<2300 




550 


-89 


Mkn 1210 




570 


-1 


<3900 




670 


+16 


NGC 4968 




270 


-26 


<2300 




400 


-31 


Circinus 




240 


-26 


< 20000 




170 


+6 


NGC 5643 




250 


-35 


<3100 




290 


+4 


NGC 1068 




640 


-27 


< 80000 




590 


-187 



Table 2. Kinematic data and Bra line decompositions from gaussian fits. For objects without BLR detection, the 
FWHM of a fit with a single gaussian is given. Total fluxes from the fits may deviate from the values in Table |l| which 
are based on a direct integration of the profile. Av describes the shift of the gaussian fitted (narrow) line centroid with 
respect to the galaxy's radial velocity as listed in Table [j]. 



the infrared that must either reflect intrinsic asymmetry 
of the NLR, or an optically extremely thick absorber that 
suppresses even infrared lines. We list in Table [2] the ve- 
locity offsets of narrow Bra and [Si IX] with respect to the 
radial velocity of Table [j] (from NED). As expected for our 
resolving power which is relatively low for detailed NLR 
profile studies, both lines are consistent with the literature 
radial velocity in most objects. The significant blueshift 
of high excitation lines in NGC 1068 has been discussed 
in detail by Lutz et al. ( |2000b[ ). In MCG-5-23-16, both 
lines are redshifted by ~80 km/s, possibly indicating an in- 
accuracy of the (optical-based) NED redshift. NGC 3081 
appears to be a case like NGC 1068, with blueshift and 
asymmetry of the coronal line persisting in the infrared. 
An interesting case deserving further study is NGC 4507, 
where a redshift of the coronal line with respect to Bra is 
suggested. 



4. Results of the Brackett a spectroscopy 

In clear cases like NGC 2992, the broad component of 
Brackett a is easily discriminated from the narrow compo- 
nent. At intermediate line widths, the situation is some- 
times ambiguous, however. It is not obvious whether a 
line originates in a true BLR of dense clouds close to 
the AGN, or whether it represents an unusually wide 
NLR profile that is seen similarly in the forbidden and 
coronal lines (see the example of NGC 1068; e.g. Lutz et 
2000bp . We hence applied the definition of a 



al. 2000a 



broad line as 'a component of Brackett a with FWHM 
around 1000 km/s or more that is not seen in the for- 
bidden or coronal lines', and used the [Si IX] line as 
reference. Table ^] summarizes our decompositions from 
gaussian fits. In cases where we do not detect a broad 
component, we quote an approximate upper limit for its 
flux determined under the assumption of a representative 
value of 3000 km/s FWHM for the BLR (e.g. Osterbrock 
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1977). One should note, however, that the spectral cov- 
erage of the data makes detection of very broad (FWHM 
10000 km/s) lines difficult even if the flux were larger. 

A number of objects deserve an individual justifica- 
tion of our decision to identify a line as originating in a 
BLR or not, also to indicate remaining ambiguities. The 
Brackett a line of Mkn 1210, for example, can be plausibly 
decomposed into a broad and a narrow component. We de- 
cided against a BLR interpretation, however, because the 
[Si IX] line shows wide wings as well, the two profiles being 
indistinguishable within the S/N limitations. This a rgues 
against the interpretation of Vcillcux et al. ( 1997 ) who 
considered Mkn 1210 a BLR detection (but see caveats in 
their appendix). Their [Fell] and Paschen (3 profiles are 
very similar to our [Si IX] and Brackett a profiles, while we 
do not find evidence for the FWHM 3000 km/s component 
possibly seen in their Brackett 7 data. 

Another complex case is NGC 5506 where a broad com- 
ponent to Pasc hen (3 has been rep orted (Blanco et al. 1990 ; 
Rixet al. |199C| ). Goodrich et al. fll994| ) and Veilleux et al. 
( |l997| ) confirm the larger width of this infrared line com- 
pared to optical lines but, on the basis of a similar [Fell] 
profile, ascribe this effect to a moderately broad but ob- 
scured component of the NLR. Our profile of Brackett a 
is relatively broad but not obviously the sum of a nar- 
row and a broad line. Comparison to [Si IX] does not help 
in this case because of its weak and uncertain detection 
in NGC 5506, making the measured FWHM highly un- 
certain. We hence compared the Brackett a profile with 
ISO spectroscopy of the [O IV] 26/im forbidden line in 
this object (Sturm et al. 2002| ) which is able to pene- 
trate large obscuring columns. We fixed the NLR FWHM 
for the ISAAC spectrum to the value of 460 km/s which 
is based on the FWHM measured from a gaussian fit to 
[OIV], corrected for the resolution difference of ISO-SWS 
and ISAAC. Adopting this NLR line width we obtain a 
residual BLR co mpone nt with FWHM similar to the one 
of Blanco et al. (|l990|) but significantly lower than the 
value reported by Rix et al. ( 1990] ), most likely reflect- 
ing line profile decomposition uncertainties of the various 
datasets. We also checked for differe nces b etween the fine- 
structure line FWH M (Sturm et al. 2002j ) and the optical 
FWHM of Veilleux Jl99l4 [l991b| ) that would be expected 
if extinction within the NLR dominates the profiles. No 
significant differences were found when convolving the op- 
tical data to the lower resolution and consistently mea- 
suring the FWHM by gaussfits. On the basis of an OI 
fluores cence line and of J band [Fell] transitions, Nagar 
et al. ( pOO^ ) argue that NGC 5506 is an obscured Narrow 
Line Seyfert 1, consistent with our decomposition which 
ascribes most of Bra to a fairly narrow BLR. 

The most uncertain case is NGC 2110 with its moder- 
ately wide Brackett a line. The nondetection of [Si IX] and 
the absence of ISO spectroscopy makes a direct compari- 
son to the forbidden lines impossible. Evidence on po ssible 
broad compon ents t o Paschen (3 is mixed (Rix et al. 199C ; 
Veilleux et al. 1997 ). We have listed this case as a narrow 



width with that of [Fell] (Veilleux et al. |1997| ). We will in- 
clude this object as uncertain when discussing the sample 
statistics below. 

Figure || summarizes the detections and limits on 
broad components of Brackett a in our sample objects. 
Relatively few broad line regions are found - only 3 to 4 
out of 13 (considering the uncertain case of NGC 2110). 
For the nondetections, limits for broad Bra components 
are of the order 1 to 2 times the narrow Bra flux. 
Assuming intrinsic presence of a BLR, as confirmed by 
spectropolarimetry for most of the objects, and an intrin- 
sic broad t o narr ow line ratio 20 equal to the median of 
the Stirpe (|l990| ) M v < -22 Seyfert Is, these limits im- 
ply an obscuration of broad Brackett a of «3 mag. For a 
Galactic extinction curve, this corresponds to an equiv- 
alent visual extinction of more than 50 magnitudes. For 
^4-4.05 A4v=0-035 (Draine 1989 ) the equivalent visu al ex- 
tinction would be 86mag, for At. 05/ Ay =0.051 (Lutz 1999 ) 
59 mag. Assuming a canonical conversion factor from vi- 
sual extinction to X-ray obscuring column of N^/Ay = 
1.79 x 10 21 cm" 2 (e.g. Predehl & Schmitt |1995| ), broad 
components to Brackett a should remain visible up to X- 
ray obscuring columns of about 10 23 cm -2 . The fact that 
all our BLR detections are below this limit is thus consis- 
tent with a Galactic value of the ratio of 4/im obscuration 
to X-ray column. 

At this point, a comparison to the factors determin- 
ing the detectability of Seyfert 2 BLRs in polarized light 
is in place. The dominant factor for BLR detectabil- 
ity in polarized light is AGN luminosity, probably both 
through the corresponding variation of host galaxy dilu- 
tion and throug h a v arying size of th e scat tering region 
(e.g. Alex ander |2001 ; Lumsden et al. 2001 ; Lumsden & 
Alexander 2001 ). The BLR detectability in polarized light 
on e mod estly depends on the X-ray column (Lumsden et 
al. 2001 ), there are BLR detections up to the highest X- 
ray columns. Since most of our objects do have detections 
of the BLR in polarized light (Table |l|), and were selected 
to be bright in [OIII], we believe that the effect of AGN 
luminosity and host dilution on the BLR detectability in 
our spectroscopy is less important than the effect of ob- 
scuration. 



5. Discussion 

Our observations have resulted in the detection of broad 
line regions in about one quarter of the Seyfert 2 sam- 
ple studied. This fraction is similar to the one reported at 
shorter near-infrared wavelengths by Veilleux et al. (1997). 
Small numbers and the difficulties discussed above to as- 
cribe moderately broad profiles for some of the objects 
to a BLR or the NLR limit a detailed intercomparison of 
these fractions. It is clear, however, that the detection rate 
at the wavelength of 4/im which samples columns up to 



line, on the basis of the similarity of the Brackett a line 



10 cm (for a Galactic extinction curve) is not much 
higher than in the shorter wavelength studies which pen- 
etrate 3 to 5 times smaller columns. This cannot reflect 
intrinsic absence of BLRs since they are detected through 
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this line horizontally, changing the adopted intrinsic ratio would move the line vertically. 



spectropolarimetry in the majority of the objects with 
4/im BLR nondetections. The majority of Seyfert 2s thus 
still has considerable BLR obscuration at 4/im. This find- 
ing is still consistent with most popular models of dust 
around AGN, though being close to constraining the ones 
that predict the most modest obscurations (e.g., Granato 
et al. 1997). It is already exceeding the modest obscura- 



Geometry effects may play a role since the samples 
are different. Maiolino et al. (2001a) mainly observed 



tions (Ay up to a few tens) predicted for most polar angles 
(9 < 90 deg) in the disk wind model of Konigl & Kartje 

dggj). 



The high 4/zm obscuration and the consistency with 
a Galactic extinction curve raise an apparent conflict 
with the abnormally low ratios of visual reddening 2?b— v 
and X-ray column in luminous AGN, as summarized by 
Maiolino et al. (2001a). Comparing reddening £b-v to- 
wards the BLR as measured from optical spectra with X- 
ray Ah, they find this ratio to be 3 to 100 times lower than 
Galactic for most of the objects in their sample. Based on 
other evidence, they suggest the ratio of visual extinction 
Ay and Ah to be lower than Galactic as well, for sev- 
eral classes of AGN. Simply lowering the dust content in 
the absorber and releasing the metals into the gas does 
not solve this discrepancy of about an order of magni- 
tude. In add ition to other effects discussed by Maiolino 
et al. ( 2001b ), a low dust-to-gas ratio should result in the 
detection of broad Brackett a lines at X-ray columns ap- 
proaching 10 24 cm -2 , but none are observed in our sample. 



Seyfert Is, Quasars and intermediate Seyferts, since op- 
tical detection of a BLR was pre-requisite for their anal- 
ysis. Here, we analyze intermediate and type 2 objects, 
that is ones viewed more 'edge-on' in the picture of a 
central torus. If, as not implausible, dust were preferen- 
tially modified or destroyed along the opening of the torus, 
then different lines of sight probe different dust properties. 
The modified dust may be found detached in the opening 
of the torus, or associated with the wa lls. Si milarly, the 
dusty wind model of Konigl & Kartje ( 1994 ) provides a 
geometry where dust-free gas is surrounded by dusty gas 
but, as noted above, predicts relatively low obscuration. 
Another geometric option is that 'normal' dust on large 
scales (100 pc or more) contributes to the obscuration, as 
well as 'modified' circumnuclear dust, the two covering 
different directions. This geometry is similar to the sce- 
nario proposed by Maiolino & Rieke (1995) for the RSA 
Seyferts, but with the difficulty that they invoke the likely 
more 'normal' large scale dust to create the obscuration 
of intermediate Seyferts. This conflicts with observations 
that intermediate Seyferts do show anomalous relations of 



optical and X-ray columns (Maiolino et al. 2001a). This 
is not consistent with 'normal' galactic dust, unless one 
postulates that there is a spatial separation between X- 
ray absorber (dust-free and close to the nucleus) and nor- 
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mal dust optical absorption on large scale (Weingartner 
& Murray pOO^ ). 

The seemingly discrepant optical and 4/im results are, 
however, also consi stent w ith the preferred interpretation 
of Maiolin o et al. ( 2001b ) for the results of Maiolino et 
al. ( 2001a ). If dust in the dense circumnuclear regions of 
AGN is domi nated by large grains, a s propo sed by Laor & 
Drainc fll993| ) and by Maiolino et al. ( j2001b| ), the net effect 
on the extinction curve will be mainly a reduction at short 
wavelengths, but much less change in the infrared. This is 
illustrated in the top panel of Fig. |j| where the Galactic 
'standard' extinction curve is compared with an extinc- 
tion curve due to a grain distribution biased in favor of 
large sizes f]. This effect would explain both the mismatch 
between optical extinction and gaseous column measured 
in the X-rays (discussed in Maiolino et al. 2001a) and the 
agreement between IR extinction and X-ray absorption 
found in this paper. 

The 'large grains' curve has a turnover at about 2- 
3/xm (at variance with the 'Galactic' curve) which is nicely 
probed by the three infrared hydrogen lines Pa/3, Br7 and 
Bra. In case of (screen) absorption, the deviation of the 
line ratios from the intrinsic case B depends significantly 
on the extinction curve, this is shown in the lower panel of 
Fig. |^. For two of the sources showing broad Bra, namely 
NGC2992 and MCG-5-23-16, there are measurements in 
the literatu re for the broad co mpone nts of Pa/3 and Br7 
(Gilli et al. |2000| , Veilleux et al. [1997D , which allow to lo- 
cate these objects in the lower panel of Fig. || The location 
of NGC 2992 is inconsistent with the Galactic extinction 
curve, but fully consistent with the 'large grains' scenario. 
For MCG-5-23-16 the location is not even consistent with 
the Galactic extinction, because of the very low Bra/Pa/3 
ratio. This might be due to variability of the BLR, since 



Bra (this paper) and Pa/3 (Veilleux et al. 1997) were not 
measured simultaneously but with a time lag of 9 years. 
This point of view is strengthened by the data of Blanco et 
al. ( 1990| ) who, another 3 years earlier, measured a 4 times 
lower broad Pa/3 flux. The caveat of variability might also 
apply to NGC 2992, though here the measurements were 
much closer in time (2 years) and the data from Gilli et 
al. (2000) have higher accuracy. Additional, simultaneous 
measurements of the broad components of the infrared hy- 
drogen lines are required to unambiguously test the 'large 
grains' scenario. 

We note that a population of large grains in the dense 
circumnuclear regions of AGN is not an ad-hoc require- 
ment, but in line with the flattened extinction curves in 
Galactic dense molecular clouds (eg. Cardelli et al. 1989) 
which are ascribed to grain coagulation. However, grain 
growth may yield the formation of complex fluffy aggre- 
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Fig. 5. Top: Changing the dust size spectrum from a stan- 
dard 'Galactic' one to one dominated by large grains may 
flatten the extinction curve in the optical and induce a 
knee in the near- infrared (see text for details). Bottom: 
Expected near-infrared broad line ratios for increasing ob- 
scuration and the two cases of standard Galactic grains 
and large grains. The adopted intrinsic ratios are for 
case B and T C =20000K, n c =10 9 c m - 3 . The two objects 
from our Bra sample for which non- simultaneous Pa/3 
and Br7 data are available are indicated. Future quasi- 
simultaneous and accurately calibrated spectroscopy is 
needed to actually execute this test. 



1 More specifically, the distribution of grain sizes has been 
modelled with n oc a - ' 3 , f3 = 2.5, a mln — 0.005/im, a max = 



Draine (1993) and by 
distri- 
bution for the Galactic dust, where /3 = 3.5, am m = 0.005/im, 
Umax = 0.25pm. See also Maiolino et al. ( 2001b| ) . 



lfim, among th ose sug gested by Laor . 
Maiolino et al. ( 2001b| ), at variance with the "standard" 



grates (e.g., Dominik & Ticlens |1997[ ) rather than the sim- 
ply larger but spherical grains, which are assumed in the 
simplified model of Fig. |5j. Finally, a dust distribution bi- 
ased for large grains may also come from the preferential 
destruction of small grains by sublimation or sputtering in 
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the circumnuclear region of AGN. Depending on the mech- 
anism producing a dust distribution that is biased toward 
large grains, the absolute extinction A^m/Na may differ. 
In case of destruction of small grains it will be below but 
similar to the Galactic value, while for coagulation at a 
fixed dust mass it will rise above the Galactic value. Fig. [| 
is consistent with modest deviations in either direction. 

The current data do not allow to discriminate between 
the 'geometrical' and the 'large grains' scenario. Both are 
plausible within the unified scenario for AGN. A better 
understanding, at least for the intermediate type Seyferts 
with moderate Ay could be gained through simultaneous 
observations of the main optical to 4/xm recombination 
lines, in order to directly trace the flattening of the ex- 
tinction curve expected in the large grain scenario. Line 
ratios formed by the broad components of Brackett a, 
Brackett 7, and Paschen (3 (Fig. |) are better suited than 
a similar diagram invoking the Balmer decrement, because 
of the reduced susceptibility to departures from case B 
(e.g. Netzer 199C and references therein). While Fig. || pro- 
poses a direct test for the presence of large grains, future 
quasi-simultaneous and accurately calibrated data will be 
needed for a conclusive result. 

6. Conclusions 

We have presented new 4/im spectroscopy of a sample 
of 12 Seyfert 2 galaxies that are well-studied in the X- 
ray, and combine these data with previous spectroscopy 
of NGC 1068. The observations are designed to probe for 
the presence of optically obscured Broad Line Regions. 
The main results are 

(i) Broad components to Brackett a are detected in 3 to 
4 of these 13 galaxies. 

(ii) The detections and limits are consistent with a 
Galactic ratio of infrared and X-ray obscuring columns. 
This result can be reconciled with the low ratios of optical 
to X-ray obscuring columns observed for several AGN if 
either the obscuring dust consists of large grains leading 
to a modified extinction curve, or if variation in dusti- 
ness or dust properties exists between directions probing 
right through the putative torus and directions closer to 
its opening. 

(iii) A survey of the coronal [Si IX] 3.94/im line shows con- 
siderable variation in its ratio to Brackett a. 

(iv) Two coronal lines of [Ca VII] and [Ca V] are detected 
for the first time in an extragalactic object, the Circinus 
galaxy. 
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